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We explore for the first time the possibilities to measure an intermediate-mass (m H = 115 - 140 GeV/c 2 ) 

' Standard-Model Higgs boson in electromagnetic proton-lead (pPb) interactions at the CERN Large Hadron 

O^l Collider (LHC) via its bb decay. Using equivalent Weizsacker- Williams photon fluxes and Higgs effective field 

(-H ■ theory for the coupling yy — » H, we obtain a leading-order cross section of the order of 0.3 pb for exclusive 

yy yy 
CZ , Higgs production in elastic (pPb — > pHPb) and semielastic (pPb — > XHPb) processes at -\fs^ = 8.8 

TeV. After applying various kinematics cuts to remove the main backgrounds (yy — > bb and misidentified 

yy — > qq events), we find that a Higgs boson with m H = 120 GeV/c 2 could be observed in the bb channel with 

a 3<r-significance integrating 300 pb -1 with an upgraded p A luminosity of 10 cm~ 2 s . We also provide for 

the first time semielastic Higgs cross sections, along with elastic tt cross sections, for electromagnetic pp, pA 

and AA collisions at the LHC . 

Q ■ PACS numbers: 14.80.Bn, 25.20.Lj 

I. INTRODUCTION 

m ' 

J> ' The standard model (SM) of particle physics predicts the existence of a scalar Higgs particle (H) to explain the 

breaking of the electroweak gauge symmetry observed in nature 11]. Direct searches for the Higgs boson at the LEP 
collider have constrained its mass above 1 14.4 GeV/c 2 at 95% confidence level (CL) 01, and global fits to precision 
electroweak data exclude nin > 154 GeV/c 2 at 95% CL Yet, it turns out that the favored intermediate-mass 
range above the LEP limit is the most difficult region for Higgs searches in pp collisions at the LHC. Indeed, 
for niH < 135 GeV/c 2 , the dominant decay mode is H — » bb with a typical cross section cr(H — » bb) » 30 pb 
at = 14 TeV |4j], which is overwhelmed by the combinatorial background from QCD b-jets production with 
cr(bb) * 500 pb. As a matter of fact, the H — > bb decay channel is now considered unaccessible 1 at the LHC 0>Hl, 
and testing the expected mass-dependent Yukawa coupling of the Higgs boson to the b-quark seems to be left open 
to study only at a next e + e~ linear collider. In the intermediate-mass range, standard Higgs searches need thus to 
resort either to rare decay modes such as H — > yy or to very stringent cuts on the final-state particles, which lead to 
jj] ■ 3-4 orders of magnitude reduction of the observed cross section. 
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In this context, the clean topologies of exclusive Higgs production in "peripheral" pp processes mediated by 
colorless exchanges - such as two gluons in a color-singlet state (Pomerons) OH or two photons |^ Q1]] - are 
attracting increasing interest lfl2i nM despite their much smaller cross sections, O(l0 4 - 10~ 5 ), compared to the 
dominant gluon-fusion or vector-boson-fusion (VBF) Higgs production channels. Exclusive events are characterized 
by wide rapidity gaps on both sides of the singly produced central system and the survival of both protons scattered 
at very low angles with respect to the beam. The final-state is thus much cleaner with just the decaying prod- 
ucts in the central detector, the signal/background is much more favorable than in parton-parton interactions, and 
the event kinematics can be constrained measuring the final protons with near-beam detectors in the LHC tunnel fl4ll . 



In this work we consider Higgs production from y-y collisions, also known as ultraperipheral collisions 
(UPCs) 11511 . with proton and ion beams at the LHC. All charges accelerated at high energies generate electromag- 
netic fields which, in the equivalent photon approximation (EPA) lllql . can be considered as (quasireal) photon 
beams 2 fl7il . A significant fraction of the pp IU3l [Till and PbPb [18] collisions at the LHC will involve for the 



* Present address at Ecole Polytechnique (Palaiseau). 

1 There are however recent developments [ 5, 6] that give some hope in the W H(bb) and ZH(bb) associated production modes. 

2 The emitted photons are almost on mass shell, with virtuality -Q 2 < 1 /R 2 , where R is the radius of the charge, i.e. Q as 0.28 GeV for protons 
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first time y-induced interactions at TeV energies with effective photon luminosities never reached before. The 
highest available photon energies are of the order of the inverse Lorentz contracted radius R of the source charge, 
Umax ~ y/R- The photon spectrum is thus harder for smaller charges, which favors proton over nuclear beams in 
the production of heavy particles. However, since the photon flux scales as the squared charge of the beam, Z 2 , 
two-photon cross sections are extremely enhanced for ion beams (Z£, = 5-10 7 for lead-lead). Particle production in 
two-photon interactions at hadronic colliders has been studied at RHIC by PHENIX IU9I1 and STAR |20l l2lll . and 
at the Tevatron by CDF JH^H. They are also part of the ALICE JH IS], ATLAS HI |H, CMS Bll and 
LHCb IT3lll32tl physics programmes with proton and/or nuclear beams. Two-photon fusion favors the production of 
spinless resonances (yy — > vector is forbidden by the Landau- Yang theorem). In this work, we propose to exploit 
the photon fluxes generated by both the proton and Pb ions at the LHC to study the possible production of the 
SM Higgs boson. Higgs photon-fusion production in proton or nucleus collisions at TeV energies is not a new 
idea (see e. g. |33tl for references). Pioneering calculations for pp IT34ll and PbPb IT35l - [37ll collisions, updated more 
recently I9lr38- l40ll . predict cross sections of a SM Higgs boson with win = 120 GeV/c 2 in the o- pp ^y y ^ H — 0.1 fb 
and crpbpb^yy^H — 10 pb ranges respectively. Unfortunately, the small value of the pp cross section and the large 
event pileup 3 in pp collisions on the one hand, and the very low design luminosities expected for PbPb, on the other, 
preclude any real measurement of the Higgs boson in those channels. 

Our proposal to study y y — > H in p A collisions presents advantages with respect to both ultraperipheral AA and 
pp collisions. First, compared to AA, one benefits from (i) beam luminosities more than 4 orders of magnitude 
larger: X. P pb ~ 10 31 cm _2 s~ 1 versus -CpbPb ~ 10 27 cm^s -1 , (ii) higher beam-beam cm. energies: -Js^ = 8.8 
TeV compared to 5.5 TeV, (iii) higher yy cm. energies (harder proton photon spectrum and smaller distance of 
approach between the centers of the radiating charges), and (v) easy removal of other photoproduction backgrounds 
characterized in the AA case by additional photon exchanges which lead to forward neutron(s) emission. The net 
result is that one can reach higher masses and yields for any centrally produced system H. The advantage with 
respect to pp UPCs is threefold: (i) a Z 2 increase in one of the photon fluxes, and the possibilities (ii) to trigger 
on and carry out the measurement with almost no pileup, and (iii) to remove most of the exclusive diffractive 
backgrounds - since the nucleus is a fragile object, Pomeron-mediated interactions in p A will, at variance with pp, 
almost always lead to the emission of a few nucleons detectable in the zero degree calorimeters. 

yy 

In the following, we present a detailed generator-level study for the exclusive pPb — > pHPb (elastic) and 

p Pb —> X HPb (semielastic) processes (Fig.Q~|i, obtained with the madgraph code supplemented with nuclear equiv- 
alent photon spectra. We compute the tree-level SM cross sections for the signal - in the Higgs effective field theory 
(HEFT) approximation - and for the expected backgrounds. We determine the expected yields in 1-year run taking 
into account the maximum attainable p Pb luminosities. We then discuss the reconstruction of the H — > bb decay for 
a Higgs boson with win = 120 GeV/c 2 , including the trigger and analysis cuts needed to minimize the backgrounds. 
After accounting for basic detector reconstruction effects (£»-jet misidentification, W?-jets invariant mass resolution), 
we determine the expected statistical significance of the measurement. Our results are promising in various fronts. 
First, they indicate that the study of the difficult //-fo-quark coupling could be accessible in this production mode 
at the LHC. Second, the observation of the yy — > H process provides an independent measurement of the Higgs-y 
coupling (likely measured previously in the traditional H —> yy discovery channel). The y y-Higgs cross section is 
generated at the one-loop level by all heavy charged particles (W and top-quark in the SM) and is thus sensitive to 
possible contributions of new charged particles with masses beyond the energy covered directly by the collider: e.g. 
via chargino and top-squark loops in supersymmetric (SUSY) extensions of the SM. 



II. CROSS SECTION EVALUATION 



The cross section for a two quasireal photon process in hadronic collisions to produce a final state H at center-of- 
mass (cm.) energy W yy (e.g. a particle H with mass win, see Fig.Q]), factorizes into the product of the elementary 



(R a 0.7 fm) and Q < 0.06 GeV for nuclei (R A ss 1.2 A [/3 fm) with mass number A > 16. 
3 Pile-up refers to the overlapping events occurring in the same bunch-crossing at high luminosity £.. For the nominal pp running at 
X = 10 34 crrT 2 s~' one expects about 25 simultaneous collisions at the LHC. 
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u, u, d, d, . 



(a)Elastic case 



(b)Semielastic case 



FIG. 1: Feynman diagrams for two-photon collisions with proton and nucleus beams producing a Higgs boson decaying into 
H — » bb. (a) Elastic production (both photons are emitted coherently and the proton and nucleus survive), (b) Semielastic 
production (the photon from the proton is emitted by a quark, the proton subsequently breaks up). 



cross section for yy — » H convoluted with the equivalent photon spectra from the two colliding beams: 
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where u)\ and a>2 are the two photon energies, and /a,b(w) are the photon fluxes at energy u) emitted by the hadrons 
A and B. The photon energies determine the cm. energy W yr = ys^r = V^i^I and the rapidity y of the produced 
system: 



W 



(2) 



For symmetric systems, the maximum effective two-photon energy W™° A occurs at y = 0, when 
( jnax _ u max x yj]j min where y = -Js^/(2mjf) is the Lorentz relativistic factor 4 and b mm * 2Ra the mini- 
mum separation between the two charges of radius Ra- Note that these "maximum" photon energies do not have 
to be interpreted as a hard cut-off but as an indication of the energy (a> > yB/b) from which the photon flux is 
exponentially suppressed. 

Table U summarizes the most relevant parameters for ultraperipheral pp, p A, and AA collisions at the LHC lfl8l 
l4lTl . It has to be noted that whereas proton-proton and nucleus-nucleus collisions are obviously part of the approved 
baseline physics progra mme of the LHC, proton-nucleus collisions are still considered as an upgrade of the heavy- 
ion programme Il4ll 14211 . Detailed studies Il42ll have nonetheless shown that is a perfectly feasible mode of operation 
of the collider, and its physics possibilities have been discussed vastly in the literature as a crucial baseline for the 
interpretation of the AA data J4U l43l Hi . At the LHC, protons and ions have to travel in the same magnetic lattice 5 
i.e. the two beams have to have the same charge-to-mass ratio Z/A. This limits the beam momentum of a given 
species to p = 7TeVZ/A for the nominal 8.3 T dipole bending field. Thus, the energy in the nucleon-nucleon 
center-of-mass is -\/T^ = 14TeV -\/(ZaZb)/(A B). In the case of AB collisions with asymmetric beam energies 
the rapidity of the cm. system shifts relative to the laboratory by Ay c m = 0.5 In [(Za B)/(Zb A)]. Thus, for pA 
collisions the rapidity shifts span a range of Ay cm = 0.35 - 0.47 for ions from oxygen to lead. For yy processes 
in ultraperipheral p A collisions, the average rapidity shift is even larger due to the harder EPA spectrum from the 
proton (e.g. {Ay) » 1.7 for the case of 120 GeV/c 2 Higgs boson, see SectionlTVli. The quoted p A luminosities have 
to be considered as unofficial but plausible values (see Section [Til Cb . The Table indicates that the "maximum" yy 
cm. energies attainable range from ^/sfy* ~ 160 GeV for PbPb to 4.5 TeV for pp (obtained taking R p = 0.7 fm). 
Two-photon fusion collisions in p Pb have JH^- 
SM Higgs boson. 



260 GeV, i.e. more than twice the most probable mass of the 



4 mpi = 0.9315 GeV/c 2 for nuclei, and m p = 0.9383 GeV/c 2 for protons. 

5 The magnetic rigidity is defined as p/Z = Br for an ion with momentum p and charge Z that would have a bending radius r in a magnetic field 
B. 
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TABLE I: Relevant parameters for photon-induced processes in A 6 collisions at the LHC: (i) beam luminosity, -Cab, (ii) nucleon- 
nucleon cm. energy, yjs m , (iii) beam energies, i^eam, (iv) Lorentz factor, y = ^s m /(2m N ), (v) effective radius (of the largest 
species), Ra, (vi) photon "cutoff energy" in the cm. frame, w raax , (vii) "maximum" photon-nucleon cm. energy, -Js™*- = 
yj2(D mia m N , (viii) "maximum" photon-photon cm. energy, Js™™, and (ix) hadronic cross section, cr,„ e / (the pp value is from l45ll . 
the p A and AA geometric cross sections are obtained with a Glauber model with cr jndNN ~ 80 mb [46]). 
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A. Equivalent photon fluxes 

In the Weizsacker- Williams approximation fl^l . the flux of equivalent photons from a relativistic particle of charge 
Z is determined from the Fourier transform of its electromagnetic field. For an extended charge with form-factor 
F(Q 2 ), such as a proton or a nucleus, the energy spectrum L/a (x) = dn y /dx, where x = a>/E is the fraction of the 
beam energy carried by the photon, can be calculated from 114711 : 



aZ 2 1 



fy/A( x ) ~ 



l/2x 2 



n 



Q - Q 



Q 4 



\F(Q 2 )\ 2 dQ 2 



(3) 



where a = 1 /137, and Q 2 is the 4-momentum transfer squared from the charge. The minimum momentum transfer 
squared, Q 2 ■ as (xm^) 2 /(l - x), is a function of x and the mass ttia of the projectile. 



For UPCs invol ving io ns it is more appropriate to calculate the spectrum of equivalent photons as a function of 
impact parameter 11481 14911 . The photon energy spectrum produced by a charge Z sweeping past a target, integrated 
on the impact parameter b from £>„„■„ to infinity, is a textbook analytical result l5(ill : 



fy/A(x) = 



aZ 2 1 



2x i *r (x;)*i(*i) - x 2 (K 2 ( Xi ) - K 2 ( Xi )) 



(4) 



where x, = xm^ b m i„, and Ko, K\ are the modified Bessel functions of the second kind of zero and first order, related 
respectively to the emission of longitudinally and transversely polarized photons. The transverse polarization 
dominates for ultrarelativistic particles (y » 1). Although this approach treats the nucleus as an idealistic hard 
sphere, the use of more realistic Woods-Saxon profiles gives effective yy luminosities only about 5% lower lf5ltl 



than the hard-sphere approximation in the range of cm. 
production of an intermediate-mass Higgs boson. 



energies W yy ~ 0.5 W"™ x dominant in the exclusive 



The figure of merit for yy processes in UPCs is £, yy = £ab d£.yy/dW y y where X.ab is the collider luminosity for 
a given A B system and d£y y /dW 7y is the photon-photon luminosity as a function of the yy cm. energy obtained 
integrating the two photon fluxes over all rapidities y, d 2 X n ,y/dW yy dy = (2 / W v v )/ v /a ( W w v /2e y ) f y / n ( W y v /2e~ y ) . 



For illustration, in Fig. [2] we show £, e ® y obtained from the parametrization of ref. 84811 of d£, yy /dW yy for ion-ion 



collisions and using the X.ab luminosities quoted 6 in Table U The curves are computed for the elastic yy fluxes. 



We note, as pointed out in (H, that the £f plot (Fig. 3) of (H], does not use correct (updated) values for £, pp at the LHC. 
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Inclusion of the semielastic fluxes would yield luminosities twice higher, as we discuss later. For comparison, we 
also plot the effective y y luminosities in e + e~ collisions at the ILC {-L e +e- = 2-10 34 cm _2 s _1 ) for -\fs = 250 GeV and 
500 GeV JH. 

The elastic p Pb two-photon luminosities at the LHC are similar to those for PbPb for low y y center of mass 
energy, W yy , and become higher for W yy > 50 GeV due to the larger p Pb beam energies. For energies of interest 
for an intermediate-mass Higgs (dotted vertical line in Fig. both yy luminosities are still almost two order of 
magnitude lower than that in proton-proton collisions. However, as we discuss in Section lTlICI there are seemingly 
no technical reasons that would prevent one to increase the instantaneous proton-nucleus luminosity by up to a factor 
(9(60) (third curve in Fig.[2]i. We remark that our study, based on madgraph, does not make direct use of the effective 
two-photon luminosities plotted in Fig. [21 although it gives similar results for the fluxes as we explain below. 




FIG. 2: Effective elastic two-photon luminosities £.AB(d£. yy I dW yy ) for pp, pPb and PbPb collisions at the LHC based on the 
beam luminosities quoted in Table[l] For pPb we show also the curve corresponding to a x60 luminosity upgrade (Sect. lTTrcl . 
The effective y y luminosities in e*e~ collisions at the ILC are also shown for -\fs = 250 GeV and 500 GeV (thin dashed lines) 115211 . 
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1. Elastic production ( p A — > pHA): 

In the case of high-energy protons, the equivalent photon spectrum can be obtained from its elastic form factors in 
the dipole approximation, F M = G 2 M and F E = (4mfc 2 E + Q 2 G 2 M )l(4m 2 p + Q 2 ) with G\ = G 2 M lpL 2 p = (1 + Q 2 /Qlr\ 
where Q 2 * 0.71 GeV 2 and p 2 p * 7.78, and reads Ha): 

fy/p(x) = ^ = I ~T ^ (X ' G »- /e o) - <P(x, QL/Ql)] , where (5) 
Xx, & - (l + «y)[-ln^ + 2^] + i ^ + c(l + |)[ln^ + ^^] (6) 

k=\ k=\ 
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withy = x 2 /(l —x). The parameters a, b and c are given by a — (1 +p 2 )/4+4m 2 /Q 2 ] « 7.16, = 1 -4ra 2 /(2o a -3.96, 
and c = (/ip - l)/£> 4 ~ 0.028. This flux 7 , with Q 2 WX ~ p 2 F , is implemented in the standard version of madgraph 
v.4 f53ll . The use of more rigorous expressions e.g. including the magnetic dipole moment of the proton fsill results 
only on small differences in the final EPA fluxes 15511 . 

In UPCs involving nuclei one could imagine, starting from Eq. @, to compute the two photons flux by simple 
multiplication and identifying b mm with the radii of the two nuclei A and B for each for photon flux. However, we 
would then miss geometrical constraints in order to have the two photons produced at the same point, outside the 
nuclei while the nuclear "halos" do not overlap. In a AB collision this implies requiring not only b\ > Ra, ^2 > Rb, 
but also \b\ - 02] > Ra + Rb- The latter condition, which excludes any overlapping configurations, prevents to derive 
a factorized formula for the y y flux in terms of the momentum fraction x\ and xo carried by each photon from their 
parents nuclei in Eq. (0|. Photon fluxes are encoded in madgraph in the same way as any parton distribution and one 
has to use factorized expressions. One option to overcome those complications, as proposed by Cahn lf56ll . would 
be to impose b\ > Ra + Rb, bo > Ra + Rb which together always exclude overlaps. However, this leads to quite 
pessimistic estimates of the joint flux 14811 since it prevents configurations where the two nucleus are very close and 
produce very energetic yy collisions. An alternative is simply not to exclude overlaps and to impose b\ > Ra and 
bi > Rg alone, i.e. to use Eq. © for both fluxes. As discussed in j48ll . in such a case, the deviation from an exact 
evaluation of the flux basically depends on the ratio of the invariant mass of the yy system over the cm. energy. For 
the production of a Higgs with mu = 120 GeV/c 2 in PbPb collisions at the LHC, one expects a deviation of a factor 
of 2. This is by far the most extreme case since the Pb nuclei are large and the available energy is reduced forcing 
the yy creation point to be very near the nuclei, including overlapping configurations. For ArAr and OO collisions, 
the corrections are already much smaller ll48Tl . A comparison of the cross sections which we obtain with factorized 
fluxes without full non-overlap condition (see later Table HI1 third column) with the results for those systems using 
the exact non-overlap condition l40ll confirms the small impact of the approximations used here. 

In the pA case, using Eq. <j4j with bj_ > Ra and Eq. Q is in fact not problematic since (i) the available energy 
is higher than in AA UPCs, (ii) the overlap between the proton and the nucleus A is reduced by the small size of 
the proton (R p as 0.7 fm) 8 , and (iii) on average, the photon radiated by the proton (for which b mm - 0.7 fm) is 
much more energetic 9 . As a consequence, the y y production point is typically far from the nucleus surface and the 
proton cannot overlap with the latter. Quantitatively, we have found that the photon momentum fraction from a Pb 
nucleus in a p Pb collisions is typically below 0.015, corresponding to bi > IRpb- Therefore, in our calculations, we 
have used the equivalent photon spectrum given by Eq. ((5]) for protons and by Eq. (O for ions with the requirement 
bmin - Ra, where the effective nuclear radii are obtained in the standard way from their mass number A via 
Ra -ro A 1 ' 3 with ro = 1.2 fm. As discussed above, in practice this condition, along with the proton flux, ensures that 
the final state is produced exclusively and outside of the colliding system p A, i.e. it avoids the hadronic overlap and 
breakup of the colliding beams. The fact that our y y cross sections agree well with other recent calculations S^l 
(see later) lends support to our approximations. 



2. Semielastic production (p A — *XHA): 

In semielastic production, Fig. \T\b), the proton does not radiate coherently. The photon flux is generated by its 
quarks and is followed by the proton breakup. As found previously f59T - l6lll for two-photon fusion processes in pp 
collisions, we expect in the pA case a comparable magnitude of the incoherent photon flux (oc A = 1) emitted by the 
quarks compared to the flux from coherent elastic emission (oc Z 2 = 1). The same is not true for UPCs involving two 
nuclei where the elastic contribution clearly dominates over the semielastic one (either from the constituent protons 

or quarks of the nuclei) due to the large Z factor in both coherent fluxes. For semielastic p A — > X HA collisions, 
one can consider the proton flux as a partonic distribution y p (x, Q 2 ), the photon being the parton, where x denotes 



7 The precise value of Q 2 m K does not actually matter since the flux is already negligible for Q 2 larger than 2 GeV 2 . 

8 Of course, the proton is not a sphere and its radius is not a well defined quantity. Effectively, electron-proton scattering fits yield a charge 
rms-radius R p = 0.89 ± 0.2 fm |57|| . and diffractive results at HERA indicate an effective proton transverse size R p = 0.65 ± 0.2 fm l58ll . 

9 Hence the larger shift in the rapidity distribution of a particle X in p A X compared to hadronic p A — > X collisions (see later). 
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the momentum fraction of the "inelastic" photon in the proton and Q~ the resolution scale at which the proton is 
probed. It can be approximated by I59ll60l,l62tl 

Q 2 ) = £ log ^ J] f -P yq (xM [«Cy. 2 2 ) + g(y, Q 2 )] (7) 

with Pyq{z) = e 2 {\ + (1 - z) 2 )/z, Qn an energy cut-off, and g(x, Q 2 ) the quark PDFs in the proton. An improved 
expression was further discussed in 16311 . In our work, we shall take advantage of madgraph features and generate 
the semielastic contribution by considering the partonic processes qy — > qH — > qbb (where q — u,u, d, d, s, s are 
all possible radiating light-quarks) convoluted with the quark PDFs in the proton 10 and the photon flux from the 
nucleus A. The contributions from all the light quarks but the u are found to amount to 1/3 of the M-quark alone, 
since they are comparatively suppressed by their smaller PDFs in the proton and/or by their lower quark-charge. 

Obviously, one has to introduce an effective mass for the quark otherwise the cross section would be logarith- 
mically divergent as Eq. (0 is if is set to 0. We take this effective mass to be m q = 300 MeV. In order to 
limit the off-shellness of the quasireal photon, we have also found reasonable to bound the maximum transverse 
momentum of the outgoing quark to p q Tmax = 5 GeV/c. This provides us at the same time with a natural value for 
the factorization scale entering the PDF for the incoming quark, Q 2 = (Pj max ) 2 - 

In the semielastic class of events, hadrons from the fragmentation of the radiating quark and from the proton 
remnants spray in the proton-direction hemisphere, while a pure y beam, unaccompanied by hadronic activity, is 
generated by the nucleus in the opposite hemisphere. To guarantee a wide enough rapidity-gap as expected in 
exclusive production, we exclude events where the jet initiated by the radiating quark q' ends up in the central 
region. All our semielastic cross sections are thus computed with the condition \y q >\ > 2.5. 



B. 77 — > H coupling (Higgs Effective field theory) 

The coupling of the scalar Higgs to photons is mostly mediated by W- and top-quark loops. The Higgs effective 
field theory (HEFT) model l65ll . where the Higgs boson couples directly to photons, can be used as an approximation 
of the Standard Model. For a not too heavy (mj < 2m,) and not very energetic {pj < 2m t ) Higgs, it is a good 
approximation to take the mass of the heavy quark in the loop to infinity. In the limit of small Higgs masses - 
below about mu = 150 GeV/c 2 which satisfies m 2 f /(4m 2 v ) < 1 - the loop induced interaction can be approximately 
described by the Lagrangian 

^yyH = -\sF tiV F^H, (8) 

where = d^Ay — dyA^ is the photon field strength tensor. Due to the Abelian nature of QED there is only one 
effective vertex g between photons and Higgs bosons. The value for the coupling constant in the HEFT model as 
implemented in madgraph is given by 

a 47/ 66 228 , 696 , 5248 4 1280 , 54528 fi 56 32 2 x 

X = 1 H T w H TZ, H TZ, H TZ, H TZ, + T° T, T7 , 

5 7tv18 v 235 1645 " 8225 w 90475 29939 1646645 w 705 987 '> 
where t, = m 2 H l{Am 2 ) and t w = m^HAm 2 ^). Higher order t, and r„, terms have been neglected. 

III. RESULTS I: CROSS SECTIONS AND RATES 

We employ the madgraph v. 4 Monte Carlo f53ll with the elastic and semielastic proton photon fluxes discussed 
in the previous Section together with the nucleus photon flux Eq. © and the HEFT model for the Higgs-photon 



We have used CTEQ6L.1 Q. 
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coupling, Eq. ©, to compute the Higgs boson cross sections in two-photon fusion processes for the systems of 
TableU The Higgs decay branching ratio to bb is obtained in madgraph with hdecay 116611 . e.g. BR(// — > bb) a 72% 
for ran = 120 GeV/c 2 . We compute also the SM cross sections for the exclusive production of bb and (possibly 
misidentified) cc and light-quark (u, d, s) pairs, which constitute the most important physical background for the 
measurement of the H — » bb channel. 



A. Signal cross sections: yy — > H 

Table HH lists the cross sections for Higgs (m# = 120 GeV/c 2 ) production in photon-photon collisions for the 
systems tabulated in Table U In Figure [3] we show our predictions for the SM Higgs production cross sections as a 
function of mu for the same systems. For m# = 120 GeV/c 2 , the cross sections span a range from 0.1 8 fb for pp up to 
18 pb for PbPb. Compared to p Pb collisions, the ratios of the Higgs cross section between the different systems are 
roughly pp : p O : p Ar : p Pb : PbPb = 1/900 : 1/50 : 1/13 : 1 : 100. It is thus apparent that the large photon flux of 
the lead ion (Z 2 ) largely compensates for the higher projectile energies of the proton or light-ion beams, as well as the 
largest photon energies attainable with the smaller species. However, when one takes into account the much smaller 
maximum luminosities at reach in the PbPb running mode, such an advantage completely disappears (see Sect. IIII CI . 

TABLE II: Production cross sections for a SM Higgs boson with m H = 120 GeV/c 2 (total and for the H — > bb decay) in elastic 

y y I . 

(A B — > AH B) ultraperipheral collisions for the LHC colliding systems listed in Table[I] For pp and p A UPCs, we also quote 
the semielastic (AB —> XH B) cross sections (with the kinematical cuts \y^\ > 2.5 and p q T < 5 GeV/c for the radiating quark). 



System 


(TeV) 


a{yy -> 
[m H = 12 
H total 


H) elastic (pb) 
GeV/c 2 ] 
H -> bb 


<x(yy -> H) 
[niH — 12C 
H total 


semielastic (pb) 
1 GeV/c 2 ] 

H^bb 


PP 


14 


0.18-10" 3 


0.13-10~ 3 


0.59-10- 3 


0.45-10~ 3 


P o 

pAi 
pPb 


9.9 
9.4 
8.8 


3.510" 3 
1.3-10- 2 
0.17 


2.5-10" 3 
9.7- 10~ 3 
0.12 


4.9- 10" 3 
1.7- 10 -2 
0.16 


3.5-10" 3 
1.410 2 
0.12 


OO 

ArAr 

PbPb 


7.0 
6.3 
5.5 


3.7- 10" 2 
0.37 
18 


2.6-10" 2 
0.26 
13 







The values we have obtained for the elastic exclusive Higgs cross sections for pp and AA (those for pA 
are published in this work for the first time) agree well with those found in the recent literature J^l H(3] ex- 
cept for the PbPb case where we overpredict the cross section by a factor two compared to l40ll due to the 
absence of the exact non-overlap condition in the convolution of fluxes as discussed previously. For pp, 
the calculation of Khoze et al. 131 makes use of the standard formula for a narrow yy resonance of spin /: 
cr(yy -> H) = 8n 2 (2J + l)T(H -> yy)/m 3 H 6(1 - m 2 /ml), which yields cr(yy -> H) - 0.12 fb, for 
inn = 120 GeV/c 2 with width T(H — > yy) ^7.9 keV/c 2 16611 . This result takes into account an effective yy 
luminosity of 1.1 ■ 10~ 3 and a gap survival factor of S 2 = 0.9, encoding the probability to produce fully exclusively 
the Higgs without any other hadronic activity from soft rescatterings between the protons. 

The signal cross sections are enhanced by a factor of two by allowing for semielastic configurations where only 
the nucleus remains intact but the proton breaks apart after the photon emission from one of its quarks (see right 
columns of Table Hill . This is the first time, as far as we can tell, that semielastic UPC Higgs cross sections appear 
in the literature. As a cross-check, we have compared the semielastic cross sections for high-mass dileptons in pp 
UPCs obtained with our madgraph prescription (see Sect. Ill A 2b with the results of ll6lll finding a good agreement. 
Note that in the pp case the cross sections are multiplied by a factor of two since there are two possibilities (one for 
each proton) to emit the photon and fragment. For p Pb, we find a Higgs cross sections for semielastic y y production, 
where the photon is radiated directly from the quarks of the proton, of cr(pPb — > yy — > X H Pb) * 0.16 fb, i.e. very 
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FIG. 3: Cross sections for the exclusive production of the SM Higgs versus its mass m H in elastic ultraperipheral nucleus-nucleus, 
proton-nucleus and proton-proton collisions at the LHC (systems listed in Tabled. 

similar to that of coherent exclusive production. We do not quote the semielastic cross-sections for nucleus-nucleus 
UPCs since in these cases the photon flux, emitted either from their constituent protons or quarks, is much smaller 
than the coherent one. 

B. Background cross sections: yy — > bb, cc, qq 

The main background to the H — > bb process is the continuum production of bb and misidentified cc and qq {q = 
u,d, s) dijets. Table Hill quotes the exclusive QQ and qq cross sections in the range of invariant masses m im , - 100 - 
140 GeV/c 2 relevant for a Higgs peak at 120 GeV/c 2 expected to be smeared by the fo-jet reconstruction resolution 
(see Section ITV D\ . Without any kinematical cut the combined bb continuum cross sections over Am - 40 GeV/c 2 
is about 25 times larger (e.g. <T bb — 8 pb for p Pb) than the Higgs cross section at 120 GeV/c 2 (see Table Hill for all 
systems. 

As we will see below, such irreducible background can be safely reduced with a few kinematics cuts. The cc 
(qq) cross sections over the same mass window are a factor of 600 (resp. 1600) larger than the Higgs signal but 
the probability of misidentifying both c-jets (resp. q{q)-]ets) as Z?-jets is only of 0.25% (resp. 2 ■ 10~ 4 ), and both 
backgrounds can be further reduced with the same selection criteria applied to remove the bb continuum. Note that, 
in agreement with our leading-order (LO) calculations for the signal, we do not consider the Next-to-Leading-Order 
(NLO) production of a heavy-quark dijet accompanied by an additional gluon radiated. This process is effectively 
eliminated by our experimental requirement of two single jets in the event. 

Table [Til] lists also the elastic top-antitop continuum inclusive cross section (for all m t j) which would be an 
interesting measurement in its own right, although as a potential H — > bb background it can be easily removed given 
the presence of two extra W decays (ff — > bb W + W~) in those events. These cross sections should be taken with a 
grain of salt since the very large yy cm. energies required (m,f ~ 340 GeV/c 2 ) are at the limit of applicability of 
our photon flux approximation. 

We obtain a semielastic cross section of continuum fr-quark dijets slightly larger (by a factor of ~30%) than the 
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TABLE III: Production cross sections for exclusive production of bb, cc and qq (q = u, d, s) with m- mv = 100 - 140 GeV/c 2 and 

_ yy - | 

for (all masses) in elastic (AB — > AqqB) ultraperipheral collisions for the LHC colliding systems listed in TableU For pp 

and p A UPCs, we also quote the semielastic (A B — » X bb B, with \y# | > 2.5 and p q T < 5 GeV/c for the photon-emitting quark) 
cross sections. 



System 


V 5 ™ 


cr(yy 


-» bb) (pb) 


a(yy -> cc) (pb) 


a(yy -> qq) (pb) 


cr(yy -> tt) (pb) 






[m bh = 100-140 GeV/c 2 ] 


[m cd = 100-140 GeV/c 2 ] 


[m^=100-140GeV/c 2 ] 


[all m ti ] 




(TeV) 


elastic 


semielastic 


elastic 


elastic 


elastic 


PP 


14 


3.4- 10~ 3 


1.1 10 2 


7.9-10-2 


0.2 


0.36- 10- 3 


pO 


9.9 


6.8-10- 2 


9.5-10 2 


1.6 


3.9 


2.7-10- 3 


pAi 


9.4 


0.27 


0.36 


6.1 


15.8 


8.1-10 3 


pPb 


8.8 


3.4 


4.5 


78. 


200. 


6.2-10-2 


OO 


7.0 


0.75 




17 


39. 


3.9- ID 2 


ArAr 


6.3 


7.6 




170 


400. 


1.010-2 


PbPb 


5.5 


420 




9.4- 10 3 


2.5- 10 4 


1.810-2 



elastic ones. We do not quote the corresponding cc and qq semielastic cross sections which, for our signal over 
background studies (Sect. UVDb . we take also as a factor of 1.3 larger than the corresponding elastic values. 

C. Event rates and p A luminosity considerations 

1. pA luminosities 

Taken at face value, the results of Fig. [3] indicate that the yy — » H cross section is maximum for lead-lead 
collisions and thus that this system should be the best suited to carry out such a measurement. Nevertheless, the 
PbPb luminosity at the LHC is severely limited mostly due to two electromagnetic processes with huge cross 
sections that affect both Pb beams ll67ll : (i) electron-positron production followed by e~ capture by one of the nucleus 

(bound-free pair production, BFPP), Pb 82+ Pb 82+ -^>Pb 82+ Pb 81+ e + , with a cross section of cr bfpp = 280 b JH, and 

(ii) Coulomb dissociation of one or both nuclei due to mutual soft photon exchange(s), 208 Pb 208 Pb-^> 208 Pb 207 Pb n 
with a cross section of o" em( j = 215 b ll69ll . Both these processes create ions with a magnetic rigidity different than 
the nominal one for 208 pfr 82 + ions, leading to beam losses and a reduction of the beam lifetime. In addition, the first 
process poses a danger of LHC magnet quenching due to the large amount of Pb M+ ions straying from the nominal 
beam orbit, impinging on and heating the superconducting dipoles. Although other technical reasons limit the 
maximum luminosities attainable with ions (see below), BFPP effects effectively reduce the maximum luminosity 
in PbPb collisions to the G>(10 27 cm _2 s- 1 ) range, i.e. seven order of magnitude below the top pp LHC luminosity, 
(5(10 34 cirrV 1 ). 

In p A collisions, the upper theoretical luminosity could be naively taken as the geometric mean of the maximum 
individual proton and ion beam luminosities, e.g. for /?Pb, £!™^> x b = yj-Cpl* ' -£™' A = 3 ■ 10 30 cm _2 s _1 . However, in 
principle the Pb-beam luminosities in p Pb could be significantly improved compared to PbPb since the y y — > e + e~ 
cross section is Z 2 = 6700 times smaller and thus the quench-limit due to BFPP could be naively raised by a <9(7000) 
factor of up to JHj££ = y/ 6700 ■ X."^ ■ X™' A = 2 ■ 10 32 cm^s-'. This is, however, only an idealistic estimate for 
several reasons. The beam-beam luminosity in a generic AB collision, X.ab, is given by the standard formula 

„ NbA N b ,B h fo y , 

£ab = F(6 C , cr , o- : ) , (9) 

An e„ B* 

where Nb is number of particles (protons or ions) per bunch in each beam, kb is the number of bunches per 
beam, fo = 11.246 kHz is the revolution frequency (given by the LHC radius), y is the Lorentz factor (the 
geometric mean of the y-factors of each beam for asymmetric systems), e„ = yjy 2 - 1 <x 2 JB* is the transverse 
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normalized emittance related to the beam size cr*, and fi* the optical function at the interaction point (IP). 
F(9 c ,<T*,cr z ) is a small reduction factor from the half-crossing angle, 9 C , and bunch length cr z , which we neglect 
in this discussion. In Table [TVl we list the beam parameters relevant for the p A running mode l4ll l42tl . We note 
that, as mentioned in the discussion of Tabled the values quoted are for now only unofficial (but plausible) estimates. 



TABLE IV: Basic beam parameters for protons and ions for p A runs at the LHC: particles/bunch Nb, number of bunches /cj, 
normalized beam emittance e„, optics p" , and associated luminosity £. AB . Possible upgraded settings (see text) may lead to a 
factor of 60 improvement of the luminosities (last column). 



System 


N t 
proton 


ion 


n 

h 

proton 


omin 
ion 


al settin 
e„ (fa 
proton 


gs 
n) 
ion 


yS* (n 
proton 


i) 

ion 


-Cab 
(crrr 2 s~') 


upgraded settings 

-Cab 
(cnrV 1 ) 


pO(9.9TeV) 


1.15-10 10 


MO 9 


2808 


592 


3.75 


1.5 


0.5 


0.5 


2.7 • 10 30 


1.6- 10 32 


pAr(9.4TeV) 


1.15-10 10 


5.5-10 8 


2808 


592 


3.75 


1.5 


0.5 


0.5 


1.5 • 10 30 


1 • 10 32 


pPb(8.8TeV) 


1.15-10 10 


7-10 7 


2808 


592 


3.75 


1.5 


0.5 


0.5 


1.5 • 10 29 


1 • 10 31 



Using these nominal beam parameters 11 and Eq. (0, we obtain the pA luminosities listed in the before-last 
column of Table [TVl The obtained default p Pb luminosity is smaller compared to the simple <9(10 30 - 10 32 crrT 2 s -1 ) 
estimates given above mainly because, conservatively, the proton intensity is reduced to 10% of its standard 
value in pp collisions l42tl . Note that the number of ions/bunch Nb are significantly lower than the proton ones, 
mainly because of space charge effects in the SPS, and intra-beam scattering limits at injection in both SPS and LHC. 

There are three potential paths to improve the p A luminosity: (i) increase the proton bunch intensity Nb to its 
standard (10-times higher) value, (ii) increment the number of Pb bunches kb (which would be possible with the 
proposed new cryogenic collimators 117011 ) by a factor of 2-3, and (iii) carry out IP upgrades which should eventually 
allow factors of two smaller p* at ATLAS and CMS. All such improvements are not unrealistic given that the 
time-scale expected for a first proton-nucleus run at the LHC is at least 4-5 years after the first pp operation. By 
then, the knowledge of the collider and the upgrades related to future LHC projects will be well advanced. The 
combined effect of such upgrades would optimistically allow one to increase the p A luminosities by a factor of 60 
(last column of Table II V I t. 



2. Higgs event rates 

The expected number of Higgs bosons expected per year in ultraperipheral p A collisions at the LHC can be 
obtained, from its production cross section (sum of elastic and inelastic channels) and the time-integrated luminosity, 
with the standard formula N = cr H ■ £ AB ■ At. The nominal LHC running time with protons (resp. ions) is 8 months 
(resp. 1 month) which, with 50% efficiency, corresponds to a run time of At » 10 7 (resp. 10 6 ) s. Using the nominal 
luminosities quoted in Table IIVI we obtain the corresponding expected Higgs events for each one of the systems 
listed in Table |V] 

With the default settings and running times, the statistics are marginal for all systems involving nuclei. For the 
nominal runs, the possibility to carry out a measurement of the Higgs boson in photon-photon collisions at the LHC 
is virtually null, except maybe for pp if one could work out a trigger - e.g. using forward proton spectrometers 1 1 411 - 
that can deal with the 20 inelastic proton-proton pileup collisions overlapping with the UPC event. A straightforward 
way to increase the expected yields by a factor of 10 would be to dedicate a full LHC year (10 7 s) to a p A run. This 
plus the upgraded luminosity settings mentioned above would readily buy one a factor of 600 increase in the expected 



For kj, and e„ one uses the smallest of the proton or nucleus values. 
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TABLE V: Expected Higgs boson (m H = 120 GeV/c 2 ) production rates per year /VHi gBS (in parenthesis those in the bb channel) in 
elastic+semielastic ultraperipheral proton-proton, proton-nucleus and elastic nucleus-nucleus collisions at the LHC, for two run 
scenarios (see text). For each system, we quote the corresponding luminosity -Cab, the running time At, and the average number 
of overlapping pileup collisions, ^/V pi | cl ,p^. 



System 




nominal runs 




upgraded p A scenario 




■Cab 


At 


^^pileup^ 


Mtiggs 


-Cab 




At 


pileup ^ 






(cirr 2 s~') 


(s) 




total (H -> bb) 


(cirr 2 s~' 


) 


(s) 




total (H -> bb) 


pp(14TeV) 


10 34 


10 7 


25 


77. (55.) 


10 34 


10 7 


25 


77. (55.) 


pO (9.9 TeV) 


2.7- 10 30 


10 6 


0.20 


0.022 (0.016) 


1.610 32 


10 7 


3.9 


13. (10.) 


pAr(9.4TeV) 


1.5-10 30 


10 6 


0.18 


0.045 (0.032) 


1-10 32 




10 7 


3.6 


30. (22.) 


pPb(8.8TeV) 


1.5-10 29 


10 6 


0.05 


0.050 (0.035) 


1-10 31 




10 7 


1 


34. (25.) 


PbPb (5.5 TeV) 


5-10 26 


10 6 


5-10~ 4 


0.009 (0.007) 


5-10 26 


10 7 


510 4 


0.15 (0.1) 



integrated luminosity (100 ptr 1 ) and, thus, of the number of Higgs counts per year (last column oflVl). Unavoidably, 
the proposed improvements of the Pb-beam intensity and /?* optics imply an enhanced probability of having various 
collisions within the same bunch crossing. The occurrence of event pileup is particularly harmful in the case of UPCs 
since it eliminates the advantages given by the clean topologies of this type of collisions. The number of overlap 
collisions can be obtained from the product of the (inelastic 12 ) reaction cross section (last column of Table HI, the 
beam luminosity (last column of Table UVb and the mean bunch distance: uVp0 eup ) = (Tab X -Lab x i^tbunch), the 
latter parameter depends on the revolution frequency of the beam and the number of bunches: (Atb unc h) = l/(/o^)- 
In the case of proton-nucleus, running with ions with the nominal bunch filling scheme, the crossing frequency is 
not well defined as not all encounters will occur at integer multiples of the 100-ns RF frequency lIToh . For p Pb with 
£ = 1.5-10 29 cm _2 s _1 , one has ^/Vpii eup ^ ~ 0.05, increasing to ~1 if we consider the luminosity upgrades discussed 
above. Clearly, the p Pb system provides the best combination of signal counting rates over pileup probability. 

IV. RESULTS II: H -> bb MEASUREMENT 

In the last section of the paper we consider in detail the possibility to measure a 120 GeV/c 2 Higgs boson 
produced in ultraperipheral pPb collisions at ^JJ^ = 8.8 TeV detected in the dominant bb decay channel 
(BR = 72%) for this mass. Our analysis is entirely based on the ATLAS/CMS detectors (trackers, calorimeters), 
needed to reconstruct the b-jets and confirm the presence of a rapidity-gap, at central rapidities (\r/ \ < 2.5). No 
additional instrumentation is needed in principle except zero-degree calorimeters to help reduce possible diffractive 
proton-nucleus interactions (see below). The generator-level rapidity (y) differential distributions for the H — » bb 
signal (histogram) and decay fo-jets (dashed histogram) are shown in Fig.|4](the semielastic distribution, not shown, 
is very similar). In our calculations we take the direction of the proton-beam coming from negative rapidities. We 
note that the yy — » //(120 GeV/c 2 ) production is peaked forward, at y — 1.7 with an r.m.s of +1 units of rapidity. 
On top of the Ay c . m . = 0.47 shift due to the asymmetric p and Pb beam energies, the proton EPA y spectrum is 
harder than the Pb one and boosts the production to even larger rapidities. The distribution of the two decay fo-jets is 
centered at the parent Higgs rapidity but it is wider (r.m.s. of + 1.4). 

The discussion presented hereafter will focus on the ATLAS and CMS experiments which feature fr-jet re- 
construction capabilities in the range needed to carry out the measurement (shaded yellow area in Fig. H|. The 
ALICE Q 

acceptance for fo-jet reconstruction is unfortunately limited to a narrow region \rj \ < 0.7 and there are 
luminosity limitations, <9(10 31 cm _2 s _1 ), linked to the latency of the Time-Projection-Chamber. We do not consider 
here either the LHCb detector J3lll . although it features good b-jet identification capabilities at forward rapidities 



We do not care here about the "harmless" elastic interactions without particle production. 



13 











a 

Q. 


pPb,^s=8.8 TeV 

= yy^H->bb 








B 










10" 3 










10- 4 










10" 5 











-10 -8 -6 -4 -2 2 4 6 8 10 

pseudorapidity r| 



FIG. 4: Rapidity-differential cross section for a Higgs boson (m H = 120 GeV/c 2 , dashed histogram) and its decay bb jets (his- 
y y - 

togram) produced in p Pb — > p if Pb with H — > fob. The (shaded) yellow area indicates the acceptance of ATLAS/CMS detectors 
for fc-jets. The dashed lines indicate the full acceptance for all other type jets. 

(77 » 2-5) and covers a fraction of the Higgs decay acceptance in p Pb collisions (provided that the proton beam 
direction points in the same direction of the apparatus). 

In order to obtain realistic estimates of the detectable number of Higgs in the bb channel, we have to consider 
the potential signal losses due to (i) the trigger efficiency (Sect. IIV Al l, (ii) the geometric acceptance of the de- 
tectors (Sect. IIV B U . (ii) the inefficiencies introduced by the finite-resolution and limitations of the experimental 
reconstruction (Sect. IIV B 2\ , and (iii) the event selection cuts aiming at removing as much as possible backgrounds 
(Sects. HVCb . We consider in detail these issues in the next subsections 

A. Triggering 

There exist detailed trigger studies for ultraperipheral AA collisions at the LHC fTtjl l30ll . The main character- 
istics of (elastic) exclusive two-photon p A events is the production of a single central particle accompanied by 
large rapidity gaps on both sides of it, and the survival of the interacting proton and nucleus. At variance with 
ultraperipheral AA collisions, one cannot make use of the Zero Degree Calorimeters (ZDC) ll7lll to tag the presence 
of a nucleus radiating a photon in ultraperipheral p A collisions. In AA UPCs, mutual Coulomb excitation of the 
incoming nuclei due to additional soft photon exchange(s) in the interaction, produce forward neutrons from the 
decay of the excited nuclear Giant-Dipole-Resonance (GDR) ll72ll . The GDR excitation probability is proportional 
to the charge squared of the incoming projectile, which is a factor of Z 2 * 6700 smaller for a proton than for a lead 
beam, and very few p Pb collisions will be accompanied by forward neutron emission. 

For most of the expected p Pb beam luminosities where no significant event pileup is expected - i.e. (Npi\ eup ^ 55 1 
proton-nucleus collision per bunch crossing - one can easily record ultraperipheral Higgs events with virtually 
zero signal loss with a level-1 (LI) trigger based e.g. on (i) two back-to-back jets, with at least one of them with 
Pt > 40 GeV/c in the central detector \ < 2.5), accompanied with an exclusivity condition given by (ii) a large 
rapidity gap Arj > 2.5 without hadronic activity along the direction of the ion emitting the exchanged photon (whose 
energy is softer than the proton one). The rates of all signal and possible backgrounds with such signatures are well 
below a few Hz, easily allocatable without any prescale within the available LI bandwidth. 

In the semielastic case, only the nucleus survives and the rapidity gap is narrower along the proton direction 
due to the presence of forward hadrons from the fragmentation of the radiating quark and proton spectator partons. 
Yet, as discussed previously, our semielastic cross sections have been computed with the requirement of no activity 
within 1 77 1 < 2.5 other than the centrally produced system. In addition, since only the proton dissociates the rapidity 
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gap in the Pb direction is unpolluted by any hadronic activity. The event topology of our H — > bb semielastic events 
is basically identical to the pure elastic production and will identically pass the LI trigger defined above. 

For the maximum luminosities, (9(10 31 ) crrr 2 s _1 , considered in the "upgraded p A scenario" one has to account 
for the possibility of one concurrent p Pb collision in the same bunch crossing as the yy interaction (see Table IVTi. 
In that case the rapidity-gap condition is not necessarily fulfilled since the exclusive event is overlapped by a normal 
(hadronic) proton-nucleus interaction. One can still use very-forward-proton detectors - like TOTEM lf73tl and 
ALFA ll74ll Roman Pots or the proposed FP420 spectrometer fl4ll - to tag at the level-2 (L2) trigger, at least, the 
surviving proton in elastic yy collisions. The nucleus cannot be tagged similarly since its momentum transfer and 
energy loss are too small to leave the nominal LHC beam env elop e. One-side tagging is a standard procedure for 
triggering on semielastic two-photon events in pp collisions 11311 . One can still count on recording UPC events 
overlapping with a hadronic collision with dedicated level-1 high-p^ Z?-jet triggers, plus a level-2 single-tagging 
of a leading proton, and separate offline the two interaction vertices (which is perfectly feasible since ATLAS and 
CMS can isolate the vertices of 25 overlapping collisions at the highest pp luminosities). In the semielastic case, 
instead of the leading-proton one can tag the jet issuing from the fragmentation of the radiating quark in the forward 
calorimeters (FCAL in ATLAS, and HF/CASTOR in CMS). All in all, such a scenario is less straightforward than 
the one considered for lower luminosities and would deserve a dedicated study beyond the scope of this exploratory 
paper. For the purpose of this study we will consider that LI and L2 triggers can be defined in p Pb collisions which 
are fully efficient with respect to the signal with counting rates in the Hz range that do not require any prescale factor. 

B. Experimental cuts 

Ideally, a complete MC simulation including parton showering and hadronization, full jet reconstruction and 
GEANT-based detector response, would give fully realistic results for our study. This is, however, beyond the scope 
of this paper. We can, however, already obtain valid estimates of the feasibility of the measurement taking into 
account the known basic detector performances, and the kinematical properties of the signal and backgrounds at the 
generator-level. 

1. b-jet acceptance: 

Full jet reconstruction in ATLAS and CMS is possible within \rj \ < 5 (or within \ij \ < 6.6 in CMS if one includes 
the CASTOR calorimeter ll75lo . Jet ^-tagging requires however tens of micrometers vertex resolutions to identify 
the secondary vertex of the decay of the leading B-meson of the jet. Such capabilities are present only within the 
\rj I < 2.5 central tracking coverage of the detectors 13 (see Fig.|4|. A realistic ATLAS/CMS cut of the type 

• both b-jets within \q \ < 2.5, 

has an acceptance of 55% (i.e. it leads to a loss of 45% of the signal) for both the elastic and semielastic components. 
The acceptance for the dominant photon-photon qq continuum background (which has a wider distribution, see later) 
is fortunately smaller (23%). One could in addition displace the vertex of the p Pb interaction point along the beam 
z-axis by up to about 0.5 m by adjusting the optics in all experiments without much loss in the luminosity |p70ll - 
Such an offset in the proton direction would represent a gain of Arj = 0.2 units of pseudorapidity for b-jets, and a 
corresponding increased geometric acceptance of 60% for H — » bb. 

2. b-jet tagging efficiency: 

Standard b-jet reconstruction in ATLAS and CMS 10, Ht 1 can be tuned to enhance either the ^-tagging efficiency 
or the purity. On the one hand, due to the large charm y y background - with cross sections e A c le 4 b = 16 times larger 
than bb - excellent ^-tagging (i.e. increased purity) is required. On the other, one needs high tagging efficiency since 



13 The combination of the CMS forward HF calorimeters and the TOTEM Tl trackers |73|| could potentially help to further extend the coverage 
for fo-jets in "particle-flow"-type analyses. 
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the bb reconstruction squares any single b-jet efficiency loss. Studying the ratio {S IB) oc e 2 b _ ta ■ o~H^bh/[ s l- m i S tag ' 
cr fg + s 2 q _ mistat , ■ cr q q] for various (correlated) (mis)tagging efficiencies from HI, leads us to choose a working 
point of 70% o-jet tagging efficiency and a mis-tag rate of 5% for c-quarks and 1.5% for light-quark and gluon jets. 
The c-quark mistagging factor is a bit optimistic but it is likely at reach with multivariate-type analyses lf7o1l after 
a few years of experience with previous pp data. Jet tagging efficiencies are slightly worst in the rapidities beyond 
\rj \ m 1.5 (for both signal and background jets), however the low particle multiplicities in two-photon events, likely 
compensate for such performance losses. In our analysis, we have therefore considered the following ^-independent 
b-jet reconstruction performances: 

• b-jet tagging efficiency: 70% for a single b-jet, 

• b-jet mistagging probabilities: 5% for a c-quark, and 1.5% for a light-quark, 

which, for our double b-jet events of interest, lead to a ~50% efficiency for the signal and a total reduction of the cc 
and qq continuum backgrounds by factors of ~400 and ~4T0 5 respectively. 

C. Background rejection 

Our ultimate goal is to have a number of Higgs events collected and a signal over background (S/B) which is 
significant enough to observe the H — > bb channel at least at the S / Vb = 3<t level. We discuss in this section 
the procedure to remove in an offline analysis any remaining heavy-quark background that could have passed the 
triggers and experimental cuts discussed above. 

1. Hadronic background: 

There may be peripheral (but still hadronic) grazing p Pb collisions with heavy-quark dijet production and wide 
77 ranges without hadronic production (above experimental thresholds) which potentially pass our trigger selection 
criteria. The very low particle multiplicity and the very small transverse momentum of the centrally produced 
system (pj a; at leading-order but smeared by the experimental resolution) expected for the very low virtuality 
of the exchanged photons in UPCs can be successfully used to distinguish them from standard peripheral hadronic 
interactions: 

• Jet multiplicity, Nj et = 2: Requiring just two b-jets within \r/ \ < 2.5 selects with 100% efficiency "clean" ex- 
clusive events and removes almost completely any hadronic interaction. This eliminates also possible genuine 
yy — > H events where an additional jet from hard gluon emission is emitted, but this is consistent with the 
tree-level cross sections presented in the paper. 

• Exclusivity: Absence of hadronic activity (above detector backgrounds) outside the reconstructed jets, within 
I77I < 2.5 for tracks and neutral particles. For neutrals, one can extend the rapidity gap condition in the Pb 
direction up to \q | < 5 (or \tj \ < 6.6 if one includes the CASTOR calorimeter in CMS), while still saving the 
semielastic signal. 

• Very low transverse momentum of the dijet system: p p ™ r < 5 GeV/c. Any photon-photon central system is 
expected to be produced almost at rest. Selecting events whose net pj is below a few GeV/c, to account for 
the experimental reconstruction of a pair of two fc-jets, eliminates any hadronic collision and still saves the 
semielastic Higgs component (with Pj <5 GeV/c, see Sect. Ill A"2l . 

The application, at the offline analysis level, of such cuts removes virtually any remaining peripheral hadronic p Pb 
which could have been recorded at the trigger-level, with zero loss of the elastic and semielastic Higgs signals. In 
addition, if needed, one can also take into account the fact that the hadronic p Pb production of bb jets is peaked at 
y « 0.47 (see Sect. [IlJ whereas the two-photon fusion events are mostly centered at y a 1.7 (Fig. |4j). 

2. Diffractive and photoproduction backgrounds: 

The experimental signatures of central diffractive (Pomeron-Pomeron or Pomeron-photon) interactions - 
exclusive central object and two rapidity gaps - are very similar to two-photon fusion processes lf77ll . Central 
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exclusive production of bb jets above pj as 40 GeV/c (IPIP — > bb) is a typical background for exclusive Higgs 
production (IPIP with cross sections of the same order as the Higgs signal itself fTill . Likewise, 

photoproduction of high-pr heavy-quarks - in photon-gluon fusion via a t- or M-channel qq pair (yg — > M>) lp78L l79tl 
or photon-Pomeron 1791 18011 processes (yZP — > £>Z?) - have event topologies (survival of the nucleus, one rapidity 
gap) similar to semielastic y y production, and cross sections in p Pb collisions of order of a few nb for b-jets above 
p T ~ 40 GeV/c, the purely diffractive contribution being -10% of it f79T,l8lll . 

Though both type of events have cross sections two to three orders of magnitude larger than our signal and they can 
potentially pass the trigger cuts and contaminate our signal sample, there are various features that separate y y-events 
from iP-induced and y-gluon,y-iP events: 

• Photon-induced interactions are less central (i.e. take place at larger impact parameters) than Pomeron-induced 
ones and, thus, the corresponding gap survival probabilities for masses (9(100 GeV/c 2 ) are much larger. The 
gap survival factor for Higgs production in pp collisions at the LHC is Sy Y = 0.9 in photon-fusion compared 
to Sjpjp = 0.03 in Pomeron-fusion fgj. In the case of proton-nucleus collisions, the situation is comparably 
much more favorable for electromagnetic Higgs production: Sy y — 0.85 versus = 8 • 10~ 4 l82ll . 

• Since the nucleus is a fragile object - the binding energy of a nucleon is just 8 Me V - even the softest Pomeron- 
mediated interactions will result in the emission of a few nucleons from the ion, detectable in the zero degree 
calorimeters. On the contrary, in purely electromagnetic p A interactions both hadrons remain intact after the 
interaction and no other forward particles are emitted. 

• The net pj of elastic y y final states is zero at LO. In the semielastic case, the typical pj of the produced Higgs 
is peaked below 1 GeV/c. In both cases, the total net pj is thus smaller than for comparable y-nucleus or 
y-proton interactions and thus is also an effective tool for separating the two classes of interactions l5Tll . 

• Due to the larger nucleus flux, heavy-quark photoproduction in p Pb is dominated by collisions of photons 
emitted by the lead nucleus with a gluon (ypbgp) or Pomeron (ypbIPp) that carry a larger fraction of the 7- 
TeV proton beam energy. The produced particles will be thus more forward boosted than our signal which 
comes from yp/, y p collisions. 

In short, inclusive high-pj- heavy-quark photon- and Pomeron-induced production will partially fill one or both 
rapidity-gaps and/or be accompanied by zero-degree neutrons, and the purely exclusive production for both 
processes has a gap-survival probability below 10~ 3 . Thus, already before any kinematics cuts, both backgrounds 
are (much) smaller than our signal. 

High- pp heavy-quark production can also take place via single (or double) resolved processes with very energetic 
photons which interact via their partonic content 118311 in collisions where the proton and the nucleus come very close 
together. We disregard this contribution in this analysis for various reasons. First, although the effective two-photon 
luminosity has indeed a large tail (Fig. [2]), the photon fluxes are exponentially decreasing and the contributions to the 
cross sections in the high-W yr region relevant for such resolved processes are very small. Second, the much more 
energetic photon flux of the proton could be potentially resolved but at the price of an interaction so close (with 
an impact parameter b likely smaller than R p + Rpb) that the nucleus would break apart after the interaction. Our 
requirement for an intact nucleus implicitly sets a limit on the probability to resolve the photon. Last but not least, 
any potential resolved photon contribution, would not only contribute to the heavy-quark background but also to the 
Higgs signal itself B] . 

3. Photon-photon continuum backgrounds: 

The only physical backgrounds to H — » bb in electromagnetic p Pb collisions are the (elastic and semielastic) 
exclusive y y — > QQ, qq processes with the cross-sections 14 quoted in Table Hill Other two-photon fusion processes, 
such as yy — > t + t~ or yy — > tt (see Table [TTTb. have final-states (e.g. particle multiplicities) different than Higgs 



The semielastic cc and qq continuum are taken to be 1.3 times the corresponding elastic cross sections as found in the bb case. 
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decay into two fr-jets. Although important (see Fig. [5]), the irreducible heavy-quark dijet backgrounds can be 
suppressed with various kinematics cuts. 



> 

2 10- 

a. 



E 

T3 

D 
T3 



10- 



10- 



pPb@^s = 8.8 TeV 



H-> bb 



y y -> bb 



Y y -> cc (x 1/400) 



100 105 110 115 120 125 130 135 140 

m inv (GeV/c) 2 (jet pair) 



FIG. 5: Invariant mass distribution of pairs of fo-jets from: (i) H — > bb events and (ii) yy — > bb and cc (misidentified with 
0.25% probability) continuum, in ultraperipheral pPb collisions at 8.8 TeV. The Higgs peak is plotted with an (arbitrary) 
Am H = 1 GeV/c 2 width. 



First, the continuum can be reduced if we require the fr-jet transverse momentum to be bigger than a significant 
fraction of the bb invariant mass. In Fig. [6] (left), we compare the pr distribution of the signal and background 
before cuts. The background is dominated by one jet with low pj (but large pi), whereas the signal peaks at 
Pt ~ mull = 60 GeV/c. Selecting events where there is at least one jet with pj = 52 -60 GeV/c removes 96% 
(97%) of the bb (cc) background while killing only 49% of the signal. Second, whereas the two Higgs decay fo-jets 
are emitted isotropically, the continuum - whose relevant Feynman diagrams have quarks propagating in the t- or 
u- channel - is peaked in the forward and backward directions (Fig. [7}. Cutting on the angle 9 between the fo-jet 
(boosted to the rest frame of the pair 15 ) and the direction of the pair (helicity frame 16 ), removes an important fraction 
of the background. The Z?-jets from the continuum are clearly peaked at |cosf?| « 1, i.e. emitted either roughly in 
the same direction as the pair or opposite to it. With a rather strict |cos#| < 0.45 cut, 82% of the continuum is 
suppressed for a 55% signal loss. 

Other cuts were tested, based e.g. on the rapidity separation \yi, - yj| between jets, without further background 
suppression power. The final set of cuts applied in our analysis is thus: 

• Transverse momentum: at least one jet with p J j' between m#/2.3 = 52 GeV/c and the kinematical limit at 
m H /2 = 60 GeV/c. 

• Acollinearity: | cos 9 \ < 0.45, where 9 is the helicity-frame angle (between the b-jet, boosted to the rest frame 
of the pair, and the direction of the pair). 

• Mass window: Invariant mass of the b-jet pairs around the Higgs mass: m pa j r = 100 - 140 GeV/c 2 (the range 
given by 2 or 3 times the width of experimental mass resolution). 



15 The b-jet is at n rads from it. 

16 One can use alternatively the Gottfried-Jackson frame, which uses the direction of the beam. The results are unchanged. 
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FIG. 6: Transverse momentum distribution of single i-jets from H — > bb and yy — > bb continuum events in ultraperipheral p Pb 
collisions at 8.8 TeV, before (pure MC-level, left) and after (right) applying all the experimental cuts discussed in the text. 
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FIG. 7: Distribution of cos(#) for fe-jets in the helicity frame from: (i) yy — > bb continuum (top histograms) and (ii) H — » bb 
signal (bottom histograms) in ultraperipheral pPb collisions at 8.8 TeV, before (pure MC-level, left) and after (right) applying all 
the experimental cuts discussed in the text. 



We note that for a given event, these selection criteria (as well as the previous acceptance and ^-tagging efficiency 
cuts) do not necessarily factorize. 



D. Signal significance 



The combined application of all experimental and background suppression cuts discussed in the previous section 
leads to a loss of about 80% of the Higgs events and an important reduction of the bb, cc and qq continuum 
backgrounds by factors of 60, 5000 and 4-10 5 respectively. In one year (10 7 s) with the upgraded 10 31 cmT 2 sT l p A 
luminosity scenario (i.e. integrating 100 pb~' of luminosity) this corresponds to a collection of N = 5 Higgs and 
N = 10 bb continuum events within 100 - 140 GeV/c 2 with small cc (N = 0.16) and negligible qq contributions. 
These values include both elastic and semielastic processes for the signal and backgrounds. The continuum 
decreases exponentially in this mass range whereas the Higgs signal peaks at run = 120 GeV/c 2 . 
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FIG. 8: Expected invariant mass distribution after 3-years of pPb running at -yJs NN = 8.8 TeV (300 pb ') for a H — > bb signal 
and residual yy — > bb continuum (histogram) after analysis cuts, fitted to a Gaussian+Exponential distribution (solid line). The 
data points are the sum of signal+background. The error bars are just statistical. The dashed line shows an exponential fit of the 
continuum alone. 

In order to determine the significance of our H — > bb signal, we generate event samples consisting of the appro- 
priate number of events after cuts in 3 years (3-10 7 s) of data-taking, i.e. 300 ptr 1 of integrated luminosity, with a 
fast Monte Carlo. We assume a bb dijet invariant mass resolution of 7 GeV/c 2 . Such a value is beyond the current 
performances for fc-jets in the range pj = 50-60 GeV/c, but may be achieved, in our underlying-event-free envi- 
ronment, with particle-flow reconstruction techniques and a good data-based knowledge of the fe-jet pj resolution 
and energy scale after a few years of LHC running. We then perform a log-likelihood fit to the pseudo-data with 
two curves: (i) one assuming that there is a signal+continuum in the region 100 < nig < 140 GeV/c 2 , and (ii) one 
assuming that there is just continuum. The significance, S , is then given by 



where Ay 2 is the difference of the x 1 of the fits for the signal and null hypotheses. We fit the peak region using a 
Gaussian fixed at m# = 120 GeV/c 2 (the Higgs mass, if any, will be already known by the time this measurement 
can be carried out) with width Aoth = 7 GeV/c 2 to match the assumed experimental mass resolution. We consider 
that the shape of the background will be known, since it can be measured with high statistics by simply removing 
the kinematics cuts applied to enhance the Higgs signal. According to our simulations, the exclusive continuum 
after cuts can be well reproduced by an exponential distribution with inverse slope [28 GeV/c 2 ] , which we also fix 
in our fit. The normalisation of the signal and background are left as the only two free fit parameters. We repeat this 
method for 500 pseudo-data sets to obtain the average significance of the fits. An example of a typical pseudo-data 
set and fits is shown in Fig. [8] The error bars correspond to "vW statistical errors. For a total integrated luminos- 
ity of 300 pb -1 , we reach a signal-to-background ratio of S H->bb/Byy^>bh * 1-5 and a statistical significance of S ~ 3. 

The main motivation for such a measurement is the unique observation of the SM H — > bb decay which seems 
otherwise unaccessible at the LHC. In addition, the observation of the yy — > H process will provide an independent 
measurement of the Higgs-y coupling, likely measured before in the H — > yy discovery channel. The yy-Higgs 
cross section is generated at the one-loop level by all heavy charged particles (W and top-quark in the SM) and it is 
thus sensitive to possible contributions of charged particles predicted in various extensions of the SM: e.g. chargino 
and top-squark loops in SUSY models, and/or charged Higgs bosons in general 2-Higgs doublet models (2HDMs). 
Last but not least, in the minimal SUSY extension of the SM - which predicts three neutral Higgs bosons: the light 




(10) 
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CP-even h, the heavy CP-even H, and the CP-odd A - the properties of the of h boson for large A masses, are similar 
to the SM Higgs boson and it could be detected in the bb decay mode as described in this work too. A dedicated 
study will be carried out to check the feasibility of such a measurement. 

V. SUMMARY 

We have presented a detailed study of the exclusive production of the SM Higgs boson in electromagnetic 

(ultraperipheral) proton-nucleus collisions at the LHC. We have evaluated the production cross sections and the 

yy yy 
corresponding yields via two-photon fusion processes in elastic (p Pb — > p H Pb) and semielastic (p Pb — > X H Pb) 

reactions, the latter being characterized by the breaking of the proton in the very forward region. Such a measure- 
ment can be used to study, on the one hand, the H - b coupling which is otherwise not accessible to measurement 
at the LHC and, on the other, the photon-photon coupling to the Higgs, i.e. provides an independent check of the 
previously measured H — » yy decay. 

First, we have computed the Higgs boson cross sections in ultraperipheral pp, p A and PbPb collisions at LHC 
energies with the madgraph Monte Carlo supplemented with equivalent photon spectra, and demonstrated that p Pb 
collisions at -\JJ^ = 8.8 TeV give the best potential for such studies when realistic reachable luminosities and event 
pile-up issues are considered. In such a case, the total cross section of a Higgs boson of ra# = 120 GeV/c 2 is about 
0.15 pb for both elastic and semielastic cases. The irreducible background due to the exclusive heavy-quark (and 
possibly misidentified light-quark) pair continuum, y y — > qq, qq has been also computed, along with the elastic it 
production cross section, for electromagnetic p Ar, p O and p Pb collisions. 

In order to determine the feasibility of the H — » bb measurement, we have proceeded to a detailed evaluation of 
the trigger setup needed, as well as the acceptances and efficiencies for the signal and continuum, and determined the 
best set of kinematical cuts needed to maximize the signal over background ratio. With reachable £>-jet experimental 
reconstruction performances, we have found that a Higgs boson with m.H =120 GeV/c 2 could be observed in the bb 
channel with a 3cr-significance integrating 300 pb~' with an upgraded p A luminosity of 10 31 cm _2 s _1 . 

To conclude, while such a study will be rather demanding in terms of luminosity, we have shown that it offers 
a unique complementary potential to the standard Higgs production mechanisms with regards to the study of its 
couplings to ^-quarks and photons at the LHC. Both measurements are key to constrain the Standard Model and any 
of its possible extensions. 
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